Correlated many-body calculation to study characteristics of Shannon information 
entropy for ultracold trapped interacting bosons. 
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We present a correlated many-body calculation to characterise Shannon information entropy of 
trapped interacting bososns. We reformulate the one-body Shannon information entropy in terms of 
one-body probability density. We approach the minimum limit of the entropy uncertainty relation 
(EUR) making N very small in our numerical work. We examine the strong correlation effect in the 
calculation of information entropy. Comparison with the mean-field result shows that the correlated 
basis function is indeed required to characterise the realistic features of information entropies. We 
also accurately calculate the point of critical instability of attractive BEC which is in close agreement 
with experimental value. Next we calculate two-body entropies in position and momentum spaces 
and introduce correaltion features to describe the strong correlation in the attractive BEC. 



I. INTRODUCTION 

In recent years, the information theoretic methods play 
an interesting and important role for the study and bet- 
ter understanding of quantum mechanical systems like 
nuclei, atomic cluster, confined atoms etc [lHg]. The syn- 
thesis of probability densities in both the position and 
momentum spaces and the associated calculation of in- 
formation entropy provide some additional relevant in- 
formation of the system. Very recently the experimental 
achievment of Bose Einstein condensation (BEC) in the 
dilute alkali atomic vapor @-E| opens a new avenue of re- 
search in this direction. The system of ultracold trapped 
atoms obey Bose Einstein statistics and a finite fraction 
of particles condense to the ground state below the crit- 
ical temperature. The accumulated atoms in the single 
state can be described as a single quantum stuff. The 
presence of external harmonic trapping potential makes 
dual possibility of measurement both in the coordinate 
and momentum spaces. The interparticle interaction and 
the confining potential strongly influence the static and 
dynamic properties of the condensate. Correlation be- 
tween the atoms in the external trap also plays a key 
role although the system is extremely dilute. 

In this paper we present the results of numerical stud- 
ies of Shannon information entropy of ultracold trapped 
interacting bosons in both the coordinate and momen- 
tum spaces. We adopt two-body correlated basis func- 
tion together with the realistic van der Waals potential 
for the description of zero-temperature condensate. From 
the condensate wave function we calculate the one-body 
density Ri(fk) and two-body density function i?2(^ij), 
where is the distance of /cth boson from the center 
of mass of the condensate and r^j is the relative separa- 
tion of the (ij) pair of bosons. The Shannon entropy in 
position space is defined as 



where p(r) is either the one-body or the pair-density in 
the coordinate space. The momentum space Shannon 
entropy is described in similar way as 
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4>(p) In cf>(p)dp 



(2) 



where 4>(p) represents the one-body or pair-momentum 
density. p(r) and ()>(p) are normalized to unity. As densi- 
ties are measured in oscillator unit (o.u.), all the entropy 
values presented here are also in o.u.. The Shannon in- 
formation entropies basically measure the uncertainty of 
the probability distribution in the respective spaces. Us- 
ing the position and momentum space entropies (S r and 
S p respectively), Bialynicki-Birula and Mycielski (BBM) 
derived a stronger version of the Heisenberg uncertainty 
relation For a three-dimensional system the total 

cntropic sum has the form 



S = S r + S p > 3(1 + In (V)) ~ 6.434. 



(3) 



p(r) In p(r)dr 



(i) 



This means that the conjugate position and momentum 
space information entropies S r and S p maintain an in- 
verse relationship with each other. It signifies that when 
a system is strongly localized in position space, the cor- 
responding information entropy and uncertainty in posi- 
tion space decreases. The corresponding momentum dis- 
tribution becomes delocalized, the information entropy 
and uncertainty in momentum space increases. It is im- 
portant to note that BBM inequality (Eq.(3)) presents a 
lower bound and the equality is maintained for a Gaus- 
sian wave function. 

It is to be noted that information theory of correlated 
bosons in the external trap has been studied using un- 
corrected mean-field theory (TlT | . The universal property 
of total entropy in different quantum systems is also pro- 
posed which takes the form S = a + bin N, a and b are 
the parameters and depend on the given system (2j. Al- 
though the choice of our system is quite similar to that 
of earlier studies, however the motivation of our present 
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work is quite different and as follows. First: as the in- 
teracting bosons in the external trap are correlated, the 
uncorrelated mean-field Gross-Pitaevskii (GP) equation 
may not reveal the correct information and characteristic 
features in the measures of Shannon entropy. In the GP 
equation, the interparticle interaction is described by a 
single parameter a s , the s-wave scattering length. How- 
ever the presence of long attractive tail in the van der 
Waal interaction will have dominating effect in the cal- 
culation of correlation function and will provide more re- 
alistic aspects. Thus our present calculation, keeping all- 
possible two-body correlation and realistic interatomic 
interaction will provide correct information and we may 
observe new features in the calculation of Shannon en- 
tropy and correlation properties of the confined bosons. 
Second: For the repulsive condensate a s is positive, with 
increase in a s or increase in number of bosons TV, the 
effective interaction parameter Na s increases. Thus for 
the repulsive condensate, with increase in Na s the system 
becomes less correlated. However for attaractive conden- 
sate a s < 0, the behavior is opposite. It is well known 
that for attractive BEC, the condensate collapses when 
TV becomes equal to some critical number N cr . Thus 
near TV = N cr , the condensate is highly correlated and 
at TV > N cr the condensate collapses. This is a very cru- 
cial point where the study of information entropy may 
provide some new characteristic features of Shannon en- 
tropy and from that we may get very accurate value of the 
stability factor which is experimentally known @. This 
is beyond the scope of study of mean-field GP equiation. 
We will also introduce correlation features to study the 
correlation properties of attractive condensate. 

This paper is organised as follows. In Sec. II we briefly 
review our theoretical approach. In Sec. Ill we present 
our calculation and results of information entropy of the 
one-body density, entropy uncertainty relation (EUR) 
and the universal property of total entropy. Sec. IV 
presents results of pair correlation and Shannon infor- 
mation entropy of two body density and statistical cor- 
relation function. 



II. METHODOLOGY 

A. Many-body calculation with correlated 
potential harmonic basis 

In this work we have employed our newly devel- 
oped correlated potential harmonic expansion methode 
(PHEM) which we have succesfully used in our earlier 
studies of different properties of condensate [l2T - flt| . Here 
we briefly describe the methode. Details can be found in 
Ref. [lMlI. 

We consider a system of A = TV + 1 identical bosons 
interacting through a two-body potential V(fij) = V(fi — 
fj) and confined in an external harmonic potential of 
frequency w. The time independent quantum many-body 



Schrodingcr equation is given by 

z—l i—1 2,/;>2 

-E]*{F 1 ,...,? A ) = 0; (4) 

where m is the mass of each boson, E the energy of 
the condensate and the trapping potential Vt ra p(ri) = 
^mui 2 rf. we can eliminate the centre of mass motion by 
using the standard Jacobi coordinates pol - [2^ defined as 

i — 2' — i i 

^ = Vi+I (?i+1 "i^^ (i = l,...N), (5) 

3=1 

The centre of mass coordinate is R = i Y^t=\ The 
relative motion of the bosons is given by 




i=l 



-E r ]^>(Ci,...,Cn) = 0, (6) 

where Vt ra p is the effective trapping potential and 

Vint((ij ■•■) Cn) is the sum of all pair- wise interactions ex- 
pressed in terms of Jacobi coordinates. Er(= E — |/t,oj) 
is the relative energy of the system. 

Hypersphcrical harmonic expansion method (HHEM) 
is an ab-initio many-body tool to solve themany-body 
Schrodingcr equation. The hyperspherical variables are 

constituted of the hyperradius r = Cf an d 

(3TV — 1) hyperangular variables which are comprised 
of 2TV spherical polar angles ("&j,<fj;j = 1,...,TV) asso- 
ciated with TV Jacobi vectors and (TV — 1) hyperangles 
(02> 03> ■•■•> 0jv) giving their relative lengths. The total 
wave function is expanded in the complete set of hyper- 
spherical harmonic (HH) basis [2l|. Therefore HHEM 
is a complete many-body approach and includes all pos- 
sible correlations. However due to large degeneracy of 
HH basis, it is highly restricted to three-body systems 
only and hence is not suitable for experimental BEC sys- 
tem containing a few thousand to few million particles. 
But the laboratory BEC is designed to be very dilute 
and the probability of three- and higher-body collision 
is negligible. This is done to avoid molecule formation. 
Therefore only two-body correlations in the many-body 
wave function and pair-interaction amongbosons are im- 
portant. This allows us to decompose the total wave 
function VP into two-body Faddccv component for inter- 
acting pair as 

A 

*=£0y(w)- (7) 

i,j>i 

It is worth to note that <f>ij is a function of two-body 
separation (r%) and the global hyperradius r. Thus the 
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effect of a two-body correlation comes through the two- 
body interaction in the expansion basis. <f>ij is symmetric 
under parity Py for bosons and satisfy the Faddeev equa- 
tion 



[T + V t 



trap 



E 



R\ <Pij 



A 

k,l>k 



(8) 



projecting on a particular PH, a set of coupled differential 
equation (CDE) for the partial wave u l K (r) is obtained 



h?_ d 2 

m dr 2 



Vtra P (r) + ^{C(C + 1) 

+4K(K + a + p + l)}-Ea\UKi(r) (15) 



where T = — ^ J^V^. is the total kinetic energy. Op- 

erating £\ - >t on both sides of Eq. (8), we get back the 
original Schrodingcr equation. In this approach, we as- 
sume that when (ij) pair interacts, the rest of the bosons 
arc inert spectators. Thus the total hyperangular mo- 
mentum quantum number as also the orbital angular 
momentum of the whole system is contributed by the 
interacting pair only. Next we expand m the sub- 
set of hyperspherical harmonics (HH) necessary for the 
expansion of V(f f ij). 



where C = I + U Kl = f K iu l K {r), a = ^ and 

p = 1 + 1/2. 

Jki is a constant and represents the overlap of the PH 
for interacting partition with the sum of PHs correspond- 
ing to all partitions (22[. The potential matrix element 
Vkk> (r) is given by 



V KK ,(r) 



(16) 



^2 V 2K+l 



m<(r) ■ (9) 



K 



Q N denotes 



the full set of hyperangles in the 3N- 
dimensional space corresponding to the (ij)th interacting 

pair and P^ff+i^iv) ^ s called the PH basis. It has an an- 
alytic expression: 



v l.rn trftf)\ 



2K+l\ 



D = 3N, 
(10) 

yo(D — 3) is the HH of order zero in the (3N — 3) dimen- 
sional space spanned by ...,(jv-i} Jacobi vectors; <j) 
is the hyperangle given by = r cos <p. For the remain- 
ing (N — 1) noninteracting bosons wc define hyperradius 

as 



Pij 



N-l 



\ K=l 



r smc 



(11) 



so that r 2 = rfj+pfj . The set of 3A— 1 quantum numbers 
of HH is now reduced to only three as for the(iV — 1) non- 
interacting pair 



h = h = •■■ — In-i — 0, 
mi = ni2 = ■■■ = tun-i = 0, 
n 2 = n 3 = ...riN-i = 0, 



(12) 
(13) 
(14) 



and 



for the interacting pair In = I, = m and 

jijv = K. Thus the 3iV dimensional Schrodingcr equation 
reduces effectively to a four dimensional equation with 
the relevant set of quantum numbers: hyperradius r, or- 
bital angular momentum quantum number I, azimuthal 
quantum number m and grand orbital quantum number 
2K + I for any N. Substituting Eq. (9) into Eq. (8) and 



B. Introduction of additional short-range 
correlation 



As the two-body interaction is represented by the s- 
wave scattering length a s only, the mean-field GP equa- 
tion completely disregards the detailed structure. A pos- 
itive value of a s represents the repulsive interaction and 
the negative a s represents attractive interaction. But in 
realistic interactions, such as the van der Waals poten- 
tial, there is always an attractive — type tail part at 

large separation and a strong repulsion at short separa- 
tion. Depending on the nature of these two parts, a s 
can either be positive or negative. In our earlier calcula- 
tions [23| we have already observed the effect of shape- 
dependent interatomic interaction in the many-body cal- 
culation. Moreover in the mean-field GP equation the 
Hamiltonian is unbound from below. So in our present 
many-body calculation, we use a realistic inter-atomic 
potential like the van der Waals potential with an attrac- 
tive — tail at large separation and a strong short range 
repulsion. The inclusion of detailed structure in the two- 
body potential with the short range repulsive core needs 
to include an additional short range correlation in the 
PH basis. This short range behavior is represented by 
a hard core of radius r c and we calculate the two-body 
wave function rj(rij) by solving the zero-energy two-body 
Schrodingcr equation 



h 2 1 d 

m r? dri- 



'' dm 



+ V(r ij ) V (r ij )=0- (17) 



This zero-energy two-body wave function ^(r^) is a good 
representation of the short range behavior of 4>ij as in the 
experimental BEC, the energy of the interacting pair is 
negligible compared with the depth of the interatomic po- 
tential. It is taken as the two-body correlation function 



4 



in the PH expansion basis. The value of r c is obtained 
from the requirement that the calculated a s has the ex- 
pected value [HI- We introduce this as a short-range 
correlation function in the expansion basis. This also 
improves largely the rate of convergence of the PH basis 
and we call it as correlated Potential Harmonic expansion 
method (CPHEM). We replace Eq. (9) by 



^Tv^mWArM^)- (18) 



K 



and the correlated PH (CPH) basis is given by 



l 2K+1 \ N )\correlated 



yl,m 



(n^Mnj), (19) 



The correlated potential matrix Vkk' (f) is now given by 



V KK ,(r) = (hfhf l )-^x 
J^{pf( z)V (r^)p^( Z )r 1 (r 



Wi{z)}dz- 
(20) 

Here P^P (£) is the Jacobi polynomial, and its norm and 
weight function are h^f and Wi(z) respectively 25 1. 

One may note that the inclusion of r\(rij ) makes the PH 
basis non-orthogonal. One may surely use the standard 
procedure for handling non-orthogonal basis. However in 
the present calculation we have checked that rj(rij) dif- 
fers from a constant value only by small amount and 

the overlap (VoK+li^N )I^2if+i(^;v )v(rkl)) is quite 



small. Thus we get back the Eq. (15) approximately 
when the correlated potential matrix is calculated by 
Eq. (20). 



III. RESULT 

A. Choice of two-body potential and calculation of 
many-body effective potential 

The interatomic potential has been chosen as the van 
der Waals potential with a hard core of radius r c , viz., 



V(nj) = oo for nj < r c and = for r,j > r c . C 6 is 

ij 

known for a specific atom and in the limit of Cq — > 0, the 
potential becomes a hard sphere and the cutoff radius 
exactly coincides with the s-wave scattering length a s . 
By utilizing the Feshbach resonance one can effectively 
tune the scattering length a s . In our choice of two-body 
potential we tune r c to reproduce the experimental scat- 
tering length. As we decrease r c , a s decreases and at a 
particular critical value of r c it passes through an infinite 
discontinuity, going from — oo to oo (l9| . For our present 
calculation we choose Rb atoms with Cq = 6.4898 x 10 -4 
o.u. [2f|. With this V(jij) we solve the zero-energy two- 
body Schrddinger equation and tune r c to obtain the a s 
correctly. We choose r c such that it corresponds to the 
zero node in the zero-energy two-body wave function for 
a s < and one node for a s > [l9|, [2(| . With these set 
of parameters we solve the coupled differential equation 



by hypcrsphcrical adiabatic approximation (HAA) [27] . 
In HAA, we assume that the hyperradial motion is slow 
compared to the hyperangular motion and the poten- 
tial matrix together with the hypercentrifugal repulsion 
is diagonalized for a fixed value of r. Thus the effec- 
tive potential for the hyperradial motion is obtained as 
a parametric function of r. We choose the lowest eigen 
potential ujq (r) [corresponding cigen column vector being 
XKo(r)} as the effective potential in which the condensate 
moves collectively. We solve the adiabatically separated 
hyperradial equation in the extreme adiabatic approxi- 
mation (EAA) 



h 2 d 2 

7T+ w oW - E R 



Co(r) = 



(21) 



subject to approximate boundary condition and obtain 
the hyper-radial wave function Co(0- F° r our numerical 
calculation we fix I = and truncate the CPH basis to 
a maximum value K = K max requiring proper conver- 
gence. Finally the many-body condensate wave function 
^ can be constructed in terms of Co( r ) and XKo{r). 

The accuracy of HAA has been tested for various sys- 
tems and it is found to be better than 1% even for long- 
range Coulomb potential. In our calculation the van der 
Waals potential has shorter range and hence HAA is ex- 
pected to be better. Moreover the confining harmonic 
potential is smooth and the hyperradial equation is com- 
pletely decoupled. This also improves the accuracy of 
HAA. Thus we can safely use HAA in our calculation 
and this reduces the numerical complications to a great 
extent. Finally the advantages of the CPHEM are as 
follows. 

i) At least for the attractive BEC, as the condensate 
becomes highly correlated, the choice of CPH basis 
(which keeps all possible two-body correlations) is 
perfectly justified for the description of correlated 
dilute BEC. 

ii) By using the HAA we basically reduced the 
3N dimensional problem into an effective one- 
dimensional problem in hyperradial space which 
provides both qualitative and quantitative picture 
of the system. As the effective quantum numbers 
are always four (for any N) we can run our many- 
body code for very high values of N. 

iii) The many-body effective potential strongly differs 
from the mean-field potential. In CP, the choice 
of contact (^-interaction in the two-body potential 
gives rise to the pathological singularity in the ef- 
fective potential. Thus the study of post-collapse 
scenario of the attractive condensate is beyond the 
scope of mean-field theory. Whereas the presence 
of short-range hard core in the van der Waals inter- 
action does not only remove the singularity, it gives 
the realistic scenario and can describe the forma- 
tion of atomic cluster after collapse. 
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B. Information entropy of one-body density 

For repulsive BEC, we consider 87 Rb atoms in the ex- 
ternal spherical harmonic trap of trap frequency 77.87 
Hz. The cut off radius r c in the van der Waals potential 
is adjusted to get the scattering length a s = 0.00433 o.u. 
= 100 a (a is the Bohr radius) which corresponds to 
the JILA experiment 0]. For attractive BEC, our choice 
is 85 Rb atoms in the JILA trap with a s = -1.832 x 10~ 4 
o.u.. From the condensate wave function \&, we calculate 
the one-body density R\{fk) [H] as 



0.16 



(22) 



The one-body density basically measures the probability 
density of finding a particle at a distance from the 
centre of mass of the condensate. The integral over the 
hypervolumc r' excludes the variable and dr' is given 
by 



dr' = r' 3N - 4 



cos^sin 3 ^" 7 , 



ir d(j>du!ijdflN-2 , (23) 



where r' is obtained from r 2 = r' 2 + 2r\ . After a lengthy 
but straight forward calculation we arrive at a close form 
given by 



/>oo /> 1 

R 1 (r k ) = V2 / T 
Jo J-i 



1 r((D-3)/2) 



_7r 3 / 2 r((Li-6)/2)J ^ r ^ 
E^o(/)x^o(r')(/^/K.)- 1 (^ /3 ^)- 1/2 <W 



KK' 



l + z 



'l-z 



D-S 



-(D-l) 

-J- 2r 2 ) dr'dz, 



where D = 3A— 3 and is the norm of the Jacobi poly- 
nomial. The integral is computed by numerical compu- 
tation using a 32-bit Gaussian quadrature. Ri(fk) basi- 
cally contains all the information of the one body density 
correlation. In Fig. 1(a) we plot the one body density 
distribution as a function of fk (in o.u.) for several num- 
ber of bososns N =50,500 and 5000. For small particle 
number the density is sharper and with increase in par- 
ticle number the density distribution gradually pushed 
out as the effective repulsion increases. Next taking the 
Fourier transformation of R\ (r/c), we obtain the one-body 
momentum distribution $i(pfc) and plot in Fig. 1(b) for 
the same set of particle numbers as in Fig. 1(a). We see 
the expected reciprocal behavior between the position 
and momentum space wave function is accordance with 
the Heisenberg's uncertainty principle. For very small N 
(~ 100), as the net effective repulsion is small, the con- 
densate wave function is close to Gaussian in both the 
spaces. However, with increase in N, when the conden- 
sate wave function in the coordinate space spreads out 



(24) 
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FIG. 1: (Color online) Plot of one body density and the cor- 
responding one-body momentum distribution for different N. 



due to increase in effective interaction Na s , the corre- 
sponding momentum space wave function squeezes ac- 
cordingly. Next to compare our many-body results with 
the mean-field we numerically solve the GP equation of 
the form 



h 2 1 i 

7T- V 2 +^mwV + g\i)(r)\ 2 ^(f) = ^(r) (25) 
Ira 2 J 



for the same set of particles as chosen in the many-body 
calculation. Here q = ivh asN is the mean-field interac- 
tion, m being the mass of the particle. The corresponding 
one-body density both in the coordinate and momentum 
spaces are presented in Fig. 2(a) and Fig. 2(b) respec- 
tively. The expected reciprocal behavior is also seen. 

Next we calculate the Shannon information entropy 
in the coordinate space using Eq. (1) and the same in 
momentum space using Eq. (2). In Table 1 we present 
the values of S r and S p and the total entropy S = S r + S p 
versus the number of particle N for 87 Rb condensate. For 
comparison we also calculate the same for the mean-field 
and present in Table 1. 
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FIG. 2: (Color online) Plot of condensate density \tp\ 2 [panel 
(a)] and the corresponding momentum distributions \4>\ 2 
[panel (b)] obtained from mean-field GP equation for different 
TV. 



TABLE I: S r , S p and S for some typical N for repulsive BEC. 



N 


GP Theory 


CPHEM 




Sp 


5 




Sp 


S 


100 


3.369 


3.087 


6.456 


3.431 


3.004 


6.434299 


500 


3.834 


2.630 


6.465 


3.782 


2.653 


6.435 


1000 


4.102 


2.398 


6.494 


4.049 


2.389 


6.437 


3000 


4.601 


1.968 


6.569 


4.544 


1.921 


6.465 


5000 


4.858 


1.750 


6.601 


4.804 


1.674 


6.478 


7000 


5.032 


1.620 


6.627 


4.987 


1.495 


6.483 



For small TV, the net effective interaction Na s is very 
small compared to the trap energy (~ hu) and the to- 
tal entropy is very close to the lower bound of the EUR. 
However with increase in particle number N, S smoothly 
increases. With increase in net repulsive interaction 
Na s , the coordinate space wave function dclocalizcs, con- 
sequently the position space entropy S r gradually in- 
creases. It signifies that the associated uncertainty in 



the coordinate space increases. The corresponding mo- 
mentum space wave function is localized, the associated 
entropy S p gradually decreases with increase in Na s . It 
indicates that the associated uncertainty in momentum 
space decreases. Although we observe same trend both 
in the many-body and mean-field results, however quan- 
titative disagreement occurs. Next to address the uni- 
versal relation for Shannon information entropy, we plot 
it as a function of In TV in Fig. 3(a) [many-body] and 
Fig. 3(b) [mean-field]. It is proposed in earlier calcula- 
tion [111 ] that for system of fermions, the universal rela- 
tion is maintained and has the form 



S = a + b\nN 



(26) 



In the mean-field results (Fig. 3b) , we retrive the straight 
line in the plot of S vs In N and the calculated values of 
a = 6.059 and b = 0.065. In the many-body calculation 
we fail to retrieve the linear relation (Eq. (26)) for the 
entire particle range. However for large particle limit 
we observe the straight increase in the total entropy as 
a function of InN. It shows that for small and finite 
particle limit the many-body calculation exhibits finite- 
size effect. Whereas for large particle limit the mean-field 
results are close to the many-body results as expected. 

As pointed earlier that for repulsive bosons, there is 
no upper limit of the number of atoms in the trap and 
the condensate is always stable for any number of atoms. 
However in case of attractive BEC, the condensate tend 
to increase its density in the centre of the trap in order 
to lower its interaction energy with increase in the atom 
number. This tendency is balanced by the zero point ki- 
netic energy which can stabilize the system. However for 
large number of bosons, the central density becomes too 
high and the kinetic energy can not balance it anymore. 
The condensate thus collapses beyond a critical number 
N cr and the corresponding stability factor is defined as 
k C r = Nc ^ as \ . In Fig. 4 we plot S r and S p as a function 
of Na s wfiere the scattering length a s is gradually tuned 
from large repulsive (when the condensate is absolutely 
stable) to attractive (when the condensate is metastable 
and then collapses at the criticality). The inverse behav- 
ior between the position and momentum space entropies 
is the consequence of the EUR. At Na s =0, we reach the 
noninteracting limit when the system behaves as ideal 
Bose gas (IBG) and the minimum uncertainty relation is 
obtained. We have checked the numerical values of S r 
and S p . They are S r = S p = 3.217; the total entropy 
S = 6.434. However for attractive BEC, with gradual 
increase in effective negative interaction, S r sharply falls 
and S p sharply increases. It indicates the condensate 
density in the coordinate space contracts near the centre 
of the trap, whereas the corresponding momentum wave 
function spreads out as expected. Our calculated value 
of the stability factor is k™ an y- bod y = 0.457, which is 
in very close agreement with the experimental value of 
k^ pt = 0.459 ± 0.012 ± 0.054 ||. In the mean-field GP 
equation stability factor is determined in the following 
way. When N < N cr , the condensate is metastable and 
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FIG. 3: (Color online) Plot of one body entropy as a function 
of In N obtained in (a) CPHEM and (b) mean-field GP theory. 

the energy functional has a local minimum (2{J. When 
N increases, the depth of local minimum decreases and 
exactly at N = N cr , the minimum vanishes and the GP 
equation has no solution. The calculated stability factor 
is k^ r p = 0.575. Thus our present many-body calcula- 
tion not only calculates the stability factor accurately, 
but also manifests the collapse of the attractive BEC as 
the simultaneous divergence in S r and S p . 

C. Information entropy of the two-body density 

The two-body Shannon information entropies in posi- 
tion and momentum spaces are defined as [30J 

S r = - J R 2 {n, ) In R 2 ( rij )dr y , (27) 

and 

S P = - J $ 2 (p i:i ) In $ 2 (p i:j )dpij (28) 

where R.2(rij) is the two-body density distribution in 
the coordinate space and $2(Pij) is the corresponding 
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Na s 

FIG. 4: (color online) Plot of one body entropy in position 
space S r and in momentum space S p against effective inter- 
action Na a . 

Fourier transformation and calculates the pair density in 
momentum space. Pair-distribution function determines 
the probability of finding (ij')th pair at a relative sep- 
aration in the coordinate space and likewise for the 
momentum space also. We calculate it as follows: 

Mm) = J J*| 2 dr" (29) 

where ^ is the many-body condensate wave function. As 
before t" excludes integration over r^. After a lengthy 
calculation we obtain closed analytic form of R2 (r^ ) [28j 

as 




Taking the Fourier transformation of i?2 ('"ij ) we get 
^2 (Pij) which gives the pair-density in momentum space. 
Now utilizing these in Eq. (27) and Eq. (28) we calculate 
Sr and Sp for various particle number keeping a s fixed 
at 0.00433 o.u. for repulsive BEC and at -1.832 x 10~ 4 
o.u. for attractive BEC. We again observe the recipro- 
cal behavior in Sr and Sp. Sr gradually increases with 
increase in Na s and the corresponding Sp gradually de- 
creases. We again observe the diverging behavior in both 
Sr and Sp near the point of critical instability [Fig. 5]. 
Now it would be interesting to compare the two-body 
information entropies with those of one-bdoy. We ob- 
serve overall similar behavior between the one-body and 
two-body quantities, however quantitative disagreement 
persists. In Fig. 4, we see that the rate of increase in S r 
is not completely balanced by the decrease in S p , which 
makes an overall slow increase in the total entropy S. 
However in Fig. 5 we observe that although Sr steeply 
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increases with Na s as in one-body, still the rate of de- 
crease of Sp is very slow. The interatomic interaction 
with a short range hard core repulsion and long range 
attractive tail plays an important role here. Specially for 
the attractive condensate, when the atoms try to form 
clusters, the strong short range repulsion comes into play. 




0.6 - 
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FIG. 5: (color online) Plot of pair-entropy in position space Sr 
and in momentum space Sp against the effective interaction 
Na e . 



D. Correlation in dilute BEC 
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FIG. 6: (color online) . Plot of the pair distribution functions 
for different N for attractive BEC. 

As pointed earlier that for attractive BEC, both the 
interatomic correlation and realistic interatomic inter- 
action plays an important role. Due to the attractive 
interaction, even in the weakly interacting gas, the ef- 
fect of correlation becomes important and we expect to 
get new physics in the study of correlation properties. 
Although the mean-field Gross-Pitaevskii (GP) equation 
with a contact (5-interaction is adequate in the discussion 
of weakly interacting Bose gas, however our correlated 
basis function is more rigorous and will provide realistic 



picture of the study of correlation properties. To char- 
acterize the effect of correlation, in Fig. 6 we plot the 
pair-distribution function i?2(?"ij) for the attractive in- 
teraction with various particle number. Pair-correlation 
vanishes as r,j — > 0, when the atoms try to form cluster 
due to strong interatomic correlation, the strong short 
range repulsion try to separate them. This realistic pic- 
ture is observed due to the use of realistic interatomic in- 
teraction. Again due to the presence of external trapping 
potential i?2 (t\j ) can not be extended beyond the size of 
the condensate. Thus i?2(?"ij) is peaked at some interme- 
diate value of Tij . Thus our results are different from the 
earlier findings of Lieb-Linger (LL) model which treats 
one-dimensional uniform Bose gas and particles interact 
via a 5-function repulsive potential [3ll . l32j . For weak 
interaction (when the number of particles are much less 
than the critical number), the correlation length is large 
which indicates weak correlation. However for particle 
number close to the critical point, the condensate be- 
comes strongly correlated, pair-correlation length sharply 
decreases and the correlation function becomes sharply 
localized. It manifests the possibility of formation of 
atomic cluster due to large two-body correlation. Next 
we calculate the healing length £ which is considered as a 
most relevant quantity to quantify the correlation of such 
highly correlated BEC near the critical point of collapse. 
It basically measures the minimum distance over which 
the order parameter can heal [H, [29|. We calculate £ 
by balancing the quantum pressure and the interaction 
energy of the condensate and plot it in Fig. 7. £ steeply 
decreases with N when the number of atoms is very close 
to the critical number. 

0.25 | , , , , 1 

0.2 
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FIG. 7: (color online) Plot of healing length £ with TV for 
attractive BEC. 

Next we calculate another useful quantity - correlation 
length L. We define L as the half-width of correlation 
function and plot it in Fig. 8. The smooth decrease in L 
with the particle number N again confirms the presence 
of high correlation in the attractive BEC. 

As the attractive condensate becomes highly correlated 
near the critical point it is quite relevant to study its 
stability by calculating the decay rates due to two-body 
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FIG. 8: (color online) Plot of correlation length L with N for 
attractive BEC. 
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FIG. 9: (color online) Plot of the loss rate Y due to two-body 
and three-body collisions for various N for attractive BEC. 

and three-body collisions. The loss rate due to two- 
body dipolar collisions and three-body recombinations 
are given by 

T = K 2 J dr\ii\ A + K 3 J dr|V| 6 (31) 

where K2 is the two-body dipolar loss rate coefficient 
and has the value (1.87 ± 0.95 ± 0.19) x 10~ 14 cm 3 /scc. 
The three-body recombination loss rate coefficient K 3 = 
(4.24t°;™ ± 0.85) x 10~ 24 cm 6 /sec ^ is the con- 

densate wave function in coordinate space and can be 
calculated as described in sec IIIB. As for attractive con- 
densate, there is a rapid increase in condensate density 
with increase in number of atoms, we may expect large 
loss rate near the critical point. For the metastable con- 
densate we calculate the total loss rate T using Eq. (31) 
as described above and plot it in Fig. 9. For small par- 
ticle number, as the condensate density is not too high, 
the probability of two or more atoms coming very close 
is relatively low. Thus the loss rate is initially small and 



then increases with the number of atoms in the trap. The 
very sharp increase in T near the critical point signifies 
that the probability of two-body and three-body colli- 
sions increases sharply which points to the possibility of 
formation of atomic cluster. 



IV. CONCLUSION 

We have reported the results of our numrerical calcu- 
lation on the Shannon information entropy of trapped 
interacting bosons and study their several characteristic 
features. We employ two-body correlated basis function 
and realistic van der Waals interaction. This is the first 
such rigorous calculation where we observe the interplay 
between the position and momentum space information 
entropies. Due to presence of external trap, the system is 
inhomogeneous and gives some additional features. We 
also numerically establish the lower bound of BBM in- 
equality which is believed to be a stronger version of 
the Heisenbcrg uncertainty principle. The observation 
of curvature in the total entropy S also differs from the 
earlier mean-field results. Our many-body results show 
that the finite size effect is prominant for small and fi- 
nite number of atoms. We also calculate the point of 
critical instability of attractive BEC. When the number 
of atoms in the trap is very close to the critical num- 
ber, the condensate becomes highly correlated and the 
mean-field GP equation is not enough to describe such 
a correlated system. Thus our calculated stability factor 
is in very close agreement with the experimental results. 
We observe the point of instability as the point where 
the information entropies in the conjugate space diverge 
simultaneously. Thus in our present calculation we not 
only calculate the instability accurately but at the same 
time we establish the point in a stronger way providing 
more information of the quantum stuff. We also calculate 
the two-body Shannon information entropies correspond- 
ing to the two-body densities in position and momentum 
spaces and discuss their similarity and differences with 
the one-body quantities. Next we specially discuss the 
correlation in dilute BEC by introducing several correla- 
tion measures like healing length £, correlation length L. 
We also calculate the two-body and three-body loss rate 
for the attractive condensate. Very sharp increase in the 
decay rates signifies that near the point ofcriticality the 
condensate becomes strongly correlated and there will 
be a possibility of formation of atomic cluster. All these 
measure signify quantitatively the presence of strong in- 
teratomic correlation in the attractive BEC which is be- 
yond the scope of studies of mean-field theory. 
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